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We present a new exactly solvable quantum problem for which the Schrodinger equation 
allows for separation of variables in oblate spheroidal coordinates. Namely, this is the quantum 
mechanical two Coulomb centers problem for the case of imaginary intercenter parameter and 
complex conjugate charges is considered. Since the potential is defined by the two-sheeted 
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arise additional possibilities in choice of boundary conditions. Detailed classification of the 
various types of boundary-value problems is given. The quasi-radial equation leads to a new 
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calculations allowing to draw conclusions about the structure of the energy spectrum are shown. 
Possible physical applications arc discussed. 
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1 Introduction 



The quantum mechanical two Coulomb centers problem {Z\&Z-i) plays a fundamental role in various 
questions of atomic physics, quantum chemistry of diatomic molecules and collision theory. The 
problem is to determine wave functions and terms of an electron moving in the field of two fixed 
charges Z\ and Z<i located at the distance R from each other. The Hamiltonian of the system, in 
atomic units (h = m e = e = 1), is given by 

I 1 ' 2 I I 1 2 I 

Due to high symmetry of (pQ), the corresponding Schrodinger equation allows for the separation 
of variables in prolate spheroidal coordinates. As early as in 1930s, Jaffe [I], Baber and Hasse 
[2] offered expressions for one-dimensional eigenfunctions in the form of series whose coefficients 
satisfy three-term recurrence relations [3]. Thereafter, based on these expansions, algorithms were 
developed to calculate the terms with relative accuracy of 1CT 12 , and the wave functions ~ 10 
[3], [I]. All of this led to considering various generalizations of the problem {Z\eZ-i) (see, for 
example, [5]). 

Application of the two Coulomb centers problem in the atomic and molecular scattering theory is 
based on the fact that the motion of the electrons and nuclei can be considered adiabatically due to 
large difference in their masses. In other words, the potential curves E(R) or quasimolecular terms 
are introduced for the colliding particles. These curves are analytic functions of the internuclear 
distance R. Transitions between two terms Ei(R) and E 2 {R) are related to their common complex 
branch point R c in the vicinity of which the energy surface looks like a corkscrew: 

AE(R) = E 2 {R) - ExiR) ~ const(i? - Rc) 1 ' 2 . 

In case of slow collisions (v <C la.e.), the branch points located near the real axis play a dominant 
role. Then, in the framework of the adiabatic approximation, it is possible to obtain a simple 
expression for the transition probability in such inelastic processes as ionization [6] . More recently, 
in the 1980s, in connection with needs of physics of thermonuclear fusion, a range of velocities 
v ps la.e. attracted interest, that is, there appeared a need to study positions of the terms singu- 
lar points in the two Coulomb centers problem, throughout the complex plane of the intercenter 
(internuclear) distance R. For the first time, such calculations were started by E. A. Solov'ev [7] 
and then continued with his coauthors [8], [9]. The main result of these studies is that the various 
types of "hidden" quasicrossings of the terms were found, and approximate expressions relating 
quasicrossing parameters to quasi-molecule characteristics and quantum numbers were obtained. 
Note that in the works [5] and [7]-[S], the consideration was conducted in prolate spheroidal coor- 
dinates. In other words, a solution of the problem (ZxeZ-i) for real R was taken as a basis for the 
generalizations, and the terms for complex R were obtained by their analytic continuation from 
the real axis. 

In this work, we consider the generalization of a different kind. Some time ago, Yu. N. Demkov 
expressed the following idea [private communication]. In going to an imaginary parameter R 
and complex-conjugate charges in the Hamiltonian (pQ), it remains Hermitian. The Schrodinger 
equation in this case will also allow a separation of variables in oblate spheroidal coordinates but 
not in prolate. Because the potential is two-sheeted, there are additional possibilities in choice of 
boundary conditions and formulation of boundary- value problems. Thus, a new exactly solvable 
quantum problem (integrable system), whose solution can not be reduced to the special case [7]-[S]. 
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The paper is organized as follows. In the next section, we briefly describe the properties of the 
potential relating to our problem. The third part deals with the separation of variables and the 
formulation of boundary-value problems. Particular attention should be paid to appearance of a 
new type of the boundary-value problems for the quasi-radial equation. Such type of the problems 
has never been considered previously, which is connected to a configuration of the equation singular 
points and configuration of the region in which the eigenfunctions are found. Description of these 
problems solution is the subject of a separate mathematical work, so in this article we restrict 
ourselves to presenting the results of numerical calculations. In the fourth part, the asymptotics of 
the eigenfunctions and terms are presented. In section 5, the results of the numerical calculations 
are given. In Conclusions we discuss possible physical applications. 

Note also that the classical analogue of our problem is known as the generalized two fixed centers 
problem and refers to celestial mechanics |10| . For the first time, the problem arose in 1961, when 
it was necessary to take into consideration the effect of non-sphericity of the Earth's gravitational 
field on satellites trajectories. 

2 The potential and its properties 

The substitution of R — > iR, Z\ — > q± + iq2, Z2 — > q\ — iq2 for the Hamiltonian (pQ), where qi, q2, R 
are real numbers, reduces it to a Hermitian operator again, inasmuch as the last two terms are 
complex-conjugate. This operator can be considered as a new Hamiltonian whose coordinate part 
we will call the potential of the generalized two Coulomb centers problem. For more clearness, we 
use the Cartesian coordinate system. Then 



It is evident that the parameters q± and q2 appear in the expression ([2]) linearly, so we represent it 
as a sum 



V = 



qi + ni 



+ 





V = V 1 (x,y,z;q 1 ,R) + V 2 {x,y, z;q 2 , R). 
Each member of the sum ([3]) has the following types of symmetry: 



(3) 



1) 



symmetry with respect to rotations around the z— axis by an arbitrary angle, 



2) 



symmetry (antisymmetry) with respect to reflection in the xy plane: 



{x,y,z) H- (x,y,-z) 



3) 



symmetry (antisymmetry) with respect to inversion: 



(x,y,z) ^ (-x,-y,-z) 



4) 



scale symmetry: 



x, y, z, R, qi, q 2 H> Ax, Xy, Xz, XR, Xq±, Xq 2 X G {M}\0. 
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The property 3) is a consequence of 1) and 2). The last property indicates that there are only two 
non-trivial parameters in the generalized quantum mechanical two Coulomb centers problem: R 
and ratio of qi/q2 or qijq\- 

Now, let us turn our attention to the fact that the expression ([2]) defines two-sheeted mapping, 
which is singular on the circle C : x 2 +y 2 = -R 2 /4, z = 0, and not at the points z± t 2 = ±R/2 as it was 
in the problem (Z\eZ-i)- Thus, the space in which the wave functions will be determined becomes 
two-sheeted, and the potential V does not already allow a simple electrostatic interpretation. 

Way out of this difficulty is the following. First, the regular branches (sheets) can be glued 
along the singular circle C to a certain analog of the Riemann surface and the potential V can be 
considered as the electrostatic potential on the extended space. Second, if we use the fact that any 
multiply connected space can be made simply connected by inserting proper barriers, and keep one 
or another branch in ([2]) fixed, we can interpret the potential V, in such space with the barrier, as 
a certain electrostatic potential. It should be emphasized that the topological considerations allow 
considerable arbitrariness in choosing form of the barrier, it is only important that its boundary 
coincides with the singular circle C. However, the requirement of the variables separation (as shown 
below) leaves only three basic variants. They are the circle C\ : x 2 + y 2 ^ -R 2 /4, z = 0, its exterior 
C2 : x 2 + y 2 ^ -R 2 /4, z = and their union C\ [j C2, that is, the entire xy— plane. 

If the branches ([2]) are placed symmetrically relative to the top and bottom sides of the barrier: 
V+ = V-, we get an analog of the simple-layer potential; when the branches are placed antisym- 
metrically: V+ = —V-, we get an analog of the double-layer potential. 

Now, let us agree on the terminology: in future, the spectral problem in the extended space 
will be referred to as two-sheeted problem, and in the ordinary space with a barrier - one-sheeted 
problem. 



3 Separation of variables and formulation of boundary- value prob- 
lems 

It is known [3] that the potentials in which the Schrodinger equation is separable in oblate spheroidal 
coordinates (£, rj, ip) must be presented in the form of: 

2 f a(Q-b( V ) cjcp) \ 

r 2 \ e + v 2 ^ (e + i)(i-v 2 )S ' 

The connection of these coordinates with Cartesian coordinates is given by the following relations: 
x= 7rV(£ 2 + 1)(1-t? 2 )cos9?, y= -V(? 2 + l)(l-?/ 2 )sin^, z = -£77, (4) 
where the variables domains D is traditionally chosen by one of the two alternative methods: 

a) £ G [0,oo), ne [-1,1], (p G [0,2tt); (5) 

b) £G (-00,00), 776 [0,1], <^g[0,2tt). (6) 

Note that the cases ([5]) and ([6]) are related ro the one-sheeted problem, since the relation (jl]) 
determines a biunivocal mapping (bijection) / : D — > M 3 . 



4 



If we proceed to consider the two-sheeted problem, the variables domain D should be chosen 
as follows: 

c) £e (-00,00), 7/ €[-1,1], ^G[0,2vr). (7) 

Then the mapping / : D — > R 3 will be single- valued, and the inverse f- 1 : R 3 ^ D - double- valued. 

Let us change the variables in ([2]) according to ((4|) and recall that the potential of the generalized 
two Coulomb centers problem can be constructed in different ways. A simple analysis shows that 
there exists a total of nine variants of such structures allowing the separation of variables, at that 
the most physically meaningful are three of them: 

v= 4(gig + g 2 ?7 sign(? ? )) 

R(e+ri 2 ) ' {) 

v _ %jj sign(Q + q 2 y) , v 

W + r/ 2 ) ' 1 J 

The variant (JSj) refers to the double-sheeted problem with the variables domain D, as well 
as to one-sheeted problem with the domain D according to ([5]) or ©. When the domain D is 
determined by ©, the impermeable barrier must be the circle C\. Then V\ is interpreted as analog 
of the simple-layer potential, and V2 - double-layer. If the domain D is determined by ([6]), the 
barrier must be the exterior of the circle C%. Then V\ is interpreted as analog of the double-layer 
potential, and V2 - simple-layer. The variants ([9j) and (jlOp are derivatives from (JSj) and refer to 
one-sheeted problems where V\ and V2 are analogs of the simple-layer potential. In the case of 0, 
the domain D is determined by ©, and in the case of (fTU|) . it is determined by ©. 

It should be noted that if we consider the generalized quantum mechanical two Coulomb centers 
problem in the half-space, the barrier must be xy— plane. The domain D is determined as follows: 

d) £e [0,00), v e [0,1], (p g [o,2tt). 



In this case, differences between (J8J), ([9]) and (I10p vanish. Let us follow the variable separation 
procedure and the formulation of boundary- value problems by the example of (|8|) , since this variant 
is basic. 

Let us represent the wave function ^j, corresponding to the term Ej(R), in the form of 

ttj = ^k qm (t V, W R) = N kqm (R)X mk (£; R)Y mq {rr, R)e im * , (11) 

where the multiindex j = {kqm} denotes the quantum number set in which k and q coincide with 
the numbers of zeros of the corresponding functions in the variables £ and 77, and the number m 
takes the values 0, ±1, ±2, .... The normalization constant Nk qm (R) is determined by the condition 



j ^tqm(^V, l P;R)^k'q' m '(tV,^R)dV = 5 kk i5 qq i5 n 
V 



where dV = ^f-(£ 2 + r] 2 )d^dr]dip - is a volume element in the oblate spheroidal coordinates. After 
substituting (|8|) and (fTT]) ) in the Schrodinger equation 

A* + 2 (E - V) * = 
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we obtain the ordinary differential system 



ar\ arj 



m 



e 2 + i 



m 



1-r? 2 



X mk {^;R) = 0, (12) 
r mg (r/;i?) = 0. (13) 



Here p 



EjR 2 



(P > 0), 



2qiR, 



-2q2R, at that p has the meaning of the 



(0 



energy parameter; a and b are the charge parameters; \ mk 
the separation constants. 



^mfc(P' a ) and = ^mq(p,b) are 



The equations (JT2J) and (113[> supplemented by the boundary conditions form the boundary value 
problems that must be solved simultaneously, and the energy spectrum can be obtained from the 
condition 



\% k (p,a) = \W(p,b). 



(14) 



The general theory of Sturm-Liouville-type one-dimensional boundary problems implies that the 
quantum numbers k, q, m remain constant for the continuous variation of the intercenter parameter 
R, and the eigenvalues A^l(p, a) or are non-degeneracy. Consequently, if the solution of 



equation (|M|) exists, it is unique. 

Let us now discuss the formulation of boundary-value problems. On the one hand, it follows 
from the most common requirements for the wave function that 



(0 



u^a-uDur^ao^)) . 



(15) 



On the other hand, in the spatial domain where the potential becomes infinite, the particle can not 
penetrate at all, that is, there must be \Pj = 0. The continuity of requires that \Pj becomes zero 
on the boundary of this domain; generally speaking, in this case the derivative of \Pj has a jump 
|llj . Thus, in the one-sheeted problem, where the circle C\ is the barrier for the wave function, we 
have: 

tf+ta =*-|ci =0. (16) 



The normal derivative must jump 



Ci 



dn. 



g(rj;qi,q 2 ,R) , 



(17) 



Ci 



where g(j]',qi,q2, R) is a smooth function of r] and parameters qi, q 2 and R. Prom the conditions 
(|15p and (|16p . we obtain the boundary conditions for quasi-radial function 



^mfc(0; R) — 0, \X mk (£; R)\ 



->• 0. 



(18) 



In the quasi-angular equation (|13p the boundary points are simultaneously singular, so to satisfy 
the condition (|15p . it should be required that the function be limited in them: 

\Y mq (±l;R)\ <oo. (19) 

If (fTHj) and (fl~9j) are satisfied, then (fl~6|) also holds, and the condition (fT7|) is satisfied automatically. 
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In the one-sheeted problem with the barrier C 2 , it is necessary to impose on the wave function 
the following condition: 

*+|c 2 =*-|cs, = 0. (20) 

For the normal derivative we obtain 



c 2 



h(rj;qi,q 2 ,R) , 



C 2 



where h(rj;qi,q2,R) is a smooth function of 77 and parameters q%, q 2 and R. From the (|20p . we 
obtain the boundary conditions for the quasi-radial and quasi-angular functions: 



\X mk (£;R)\ 



£— >±oo 



> 0, Y mq (0; R) = 0, \Y mq {+l- R)\ < 00 



(21) 



In the two-sheeted problem no barriers appear, but there is a circle C on which the potential is 
singular, therefore, the following condition must be satisfied: 

*|cr = ¥(o,o,¥>;fl) = 0, (22) 

which is in contradiction with the expression (|lip . Indeed, if the expressions (jlip and (|22|) are 
valid, we have the alternative: 



1. X mk {0; i?) = , or Y mq {0; R) = 0; 

2. X mk (0; R) = and Y mq (0; R) = simultaneously, 

which, generally speaking, does not follow from any physical considerations. The contradiction 
arises particularly acute when we consider the ground state. 

It is apparent that both cases are not satisfied. The way out of the situation is to separate 
factors (factorization) 

*i = N kqm (R)(f + ri 2 )X mk {i- R)Y mq (rj; R)e im * . (23) 

Here, the overline means that we are dealing with one-dimensional functions of the two-sheeted 
problem. Although the functions satisfy (|12p and (|13p. but they essentially different from the 
corresponding functions of the one-sheeted problem. 

The boundary conditions in the two-sheeted problem are formulated as follows: 

\X m k(Z;R)\ —— ->0, \Y mq (±l;R)\ <oo, (24) 

because now the boundary points are singular. 

Let us now discuss the specific of the boundary-value problems related to the quasi-radial 
equation. In the problem {Z\eZ-i) the corresponding equation was considered on the interval 
[l,oo), that is, between the singular points. Then, using the Jaffe transformation [1] for the 
independent variable x = (£ — l)/(£ + 1), the interval is transferred to the unit segment, and the 
additional singular point -1 goes to —00 and does not effect the convergence of series representing 
the eigenfunction. In our problem, the equation (|12p on the complex £— plane has three singular 
points: two finite regular point £1 = +1, £2 = — 1 an d one irregular at infinity. 

It is clear that the points £1 and £2 are "on different sides" of the real axis, and so any attempt 
to present X mk (£,;R) in the form of series faces the problem of the circle of convergence. Use 
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of standard techniques, such as the transformation of Jaffe [I] or Jaffe-Lay [12], are no longer 
possible here, and therefore it is necessary to look for a more general method including the standard 
methods as special cases. Note that the problem of the circle of convergence will always appear 
in the boundary- value problems after separation of variables in the Schrodinger equation in oblate 
spheroidal coordinates (it is typical for these coordinates). 

In conclusion of this paragraph we note that the separation of the factors in (123f) . without 
changing the structure of the equations (fT2j) and (fT3j) . can affect the convergence of the series for 
the one-dimensional functions. 



4 Asymptotic behavior of the eigenfunctions and terms as R — > 0. 

In the context of possible practical applications of the model under consideration, the pattern of 
the terms is of the most interest. In order to correctly identify the individual terms it is necessary 
to know the initial approximation or asymptotic of Ej(R) as R — > 0. In this limit, the differ- 
ence between the one-sheeted and two-sheeted problems disappears and, in addition, the spherical 
symmetry arises, therefore it is natural to go to spherical coordinates (r, tp). 

Let us expand the potential ([2]) in powers of the small parameter (R/r) using the formula (see, 
for example, [33]): 

— - = y^t k P k (u) , \t\ < min\u ± \/u 2 - 1| , 
Vl-2tu + t 2 ~ 

where Pk(u) is Legendre polynomials. As a result of the transformations we obtain the expression 

^ = ^E(-ir(^) 2n ^(cos^) + ^f;(-ir i (^) 2n+1 p 2n+1 (cos^), (25 ) 

where odd multipole moments associated with the parameter q\, and the even multipoles with q<i- 
It is clear that only one term 2q\/r in (|25p does not vanish as R — > 0. This term is naturally 
chosen as the unperturbed potential while the perturbation will be the first term of the second 
sum —q2R/r 2 cos , d. Then the terms Ej(R) must continuously transit to N 2 — fold degenerate energy 
levels of the hydrogen atom with a charge 2q\ : 

E ^T^ Elflm = ~W 1 (26) 

where i = {Nlm}— is a set of spherical quantum numbers which are related to the numbers of zeros 
of the one-dimensional eigenfunctions by the relations: 

N = k + q + m + l, I = q + m. 

The radial X m k(C; R) and angular Y mq (rj;R) parts of the eigenfunction are reduced to into the 
radial Rxi(r) and angular Yj m (#,</?) parts of one-center problem. 

VkgmfoWiO) = N kqm (0)X mk (Z;0)Y mq ( m 0)e im(fi = NR Nl (r)Y lm (tf,ip) , (27) 

that will be the correct function of the zero approximation. Arguments similar to those in [3], give 
the following expression for the energy in the first non-vanishing order: 

F (nam- 2<1 " 8q 2 (q 2 + q 2 )[l(l + 1) - 3m 2 ]R 2 2 

E N i m (qi, Q2,R)-- W - m(l + m _ 1)(2Z + 1)(2Z + 3) + 0((R) ). (28) 
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Note that formula (I28p reduces to Baber and Hasse's result [2] by means of the substitution q\ — > 
[Z\ + Z 2 )/2, q 2 — > i{Z\ — Z 2 )/2. The calculation of the energy in the following order of the 
perturbation theory becomes quite cumbersome and, in addition, there one can expect logarithmic 
corrections similar to those obtained in |14j for the problem (ZieZ 2 ). 



5 The results of numerical calculations 

The boundary-value problems, which are equations fjl2[> and fjl3f) supplemented by the boundary 
conditions (fl"8|) . (fT9|) and (|2"T]). have been solved numerically. As noted above, the condition (fT4"|) 
for the fixed parameters a, b, and j = {kqm} has a unique solution 

p*{2q 1 R,-2q 2 R) = ^yf=2E~. (29) 

Solving it for Ej, we find discrete spectrum. For the classification of the terms, in addition to 
j = {kqm}, we will also use the set of spherical quantum numbers i = {Nlm}. Let us keep the 
traditional spectroscopic notation, when the numbers 

1 = 0, 1, 2, 3, 4, 5, . . . and m = 0, 1, 2, 3, . . . 

correspond to the letter rows 

L = s, p, d, f, g, h, . . . and M = o, it, 5, <f>, ■ ■ ■ 

Therefore, in order to make any conclusions about the structure of the energy spectrum in the 
generalized quantum mechanical two Coulomb centers problem, at least on a qualitative level, it is 
necessary to find a meaningful collection of the special cases of Ej{R) which covers the transition 
and asymptotic (as R — > , R — )■ oo) domains. 

Calculations with different combinations of q± = 1 , ... 10 and q 2 = 1 , . . . 10, show that such 
collection is formed from the first ten Ej(R) at R £ [0;20]. It can be explained by the fact that 
there is a scaling (see property 4) of the potential ([2])) in our problem. Thus, in spite of the fact 
that the total of about two thousand curves have been calculated, the structure of the spectrum 
was determined in a very wide domain of variation of q\, q 2 and R. 

Analysis of the curves shows that the pattern of the terms varies considerably for q\ < q 2 and 
q± > q 2 . This statement becomes especially obvious when selecting two systems differing from each 
other by the permutation q\ -H- q 2 , for example, q± = 1, q 2 = 3 (see Fig. [[]) and q\ = 3, q 2 = 1 (see 
Fig. ED. 

First of all, pay attention to the fact that the terms in Fig. [1] have distinct local minima at 
finite R, whereas the terms in Fig. [2] change rather aslope. The availability of the minima of Ej{R) 
at q\ < q 2 indicates the stable states of a moving charged particle in such systems. 

The case of q\ = requires a separate consideration. It is well-known [15] that in the spherically 
symmetric field Cr~ 2 with C ^ —1/4 the energy levels exponentially condense to the continuous 
spectrum border. One can see a similar effect in the Fig. [3j We stress that the continuous spectrum 
starts at finite values of R. 

Note that in all three figures the dashed lines depict the border of continuous spectrum. 

The question of the configuration interaction of the terms or their quasicrossings requires special 
consideration. In the works of E. A. Solov'ev [7]-[9], the branch points were found near the 
imaginary axis R, so we can expect quasicrossings in the spectrum of our problem. 
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6 Conclusion 



We have considered a new exactly solvable quantum problem. In distinction to the {Z\eZ-i) problem, 
our problem allows for various choices of boundary conditions. It gives us additional opportunities 
in setting the boundary problems. We have shown that the boundary problems associated with 
the quasi-radial equation belong to a new type. The peculiarity of these problems is related to the 
fact that the singular points of the differential equation are on different sides of the region in which 
the eigenfunctions are defined. Analytic representations of such eigenfunctions are not known yet. 
Even in the famous paper [16] where an analogous problem for spheroidal equations has arised, 
the boundary problems were not duly considered and the relevant discussions were largely avoided. 
Therefore, our paper can be considered as an incentive for the mathematical physics experts to 
analyse the problems of this type. 

We have shown for the first time that quantum problems can be considered on Riemannian 
surfaces. Probably, this idea could give some insights into the puzzle of confinement of quarks. At 
least, one could consider models of [T7], |18j and |19| type on Riemannian surfaces. 

Numerical calculations have shown that for q± < q2 the terms Ej(R) do have minima at finite 
values of R. It implies that such systems possess stable states. We believe that the Demkov problem 
may be viewed as a model for higher excited electron states in the fields of various ring molecules 
(e.g., aromatic hydrocarbons). In particular, the asymptotic formulas (]27p and (I28p may be of use 
for description of higher excited states of benzene [20] which have been observed experimentally. 
Exact quantum mechanical computation of these states is very challenging from the technical side. 
Our work offers a possibility of simple description not only for benzene but for a good number of 
ring molecules. 
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Figure 3: Terms of system q\ = 0, q<i = 1. 
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